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We discuss a possible test of neutrino oscillation hypothesis by proposing the combined analysis of 
high energy atmospheric neutrino induced muon events that have been detected around horizontal 
direction in the Kolar Gold Field (KGF) underground site and below the horizontal direction by 
many large detectors such as Super-Kamiokande and MACRO. Up/down asymmetry obtained using 
contained events recorded by detectors at Kamioka site probes low energy region of atmospheric 
neutrino whereas, the suggested method probes high energy neutrinos. It mainly depends on the 
observations and it is free of uncertainties in neutrino flux, interaction cross section etc. In this paper 
we demonstrate that the method is sensitive to a region of oscillation parameter space that explains 
all the features of atmospheric neutrino data in the Super-Kamiokande detector; the limiting factor 
being the statistical strength of the KGF observations. This method provides the only way to study 
the up/down asymmetry beyond Multi-GeV region which is yet to be measured experimentally. 



I. MOTIVATION 

Evidence for oscillation of atmospheric neutrino has 
been seen in Super-Kamiokande detector (SK) (lj. For a 
successful interpretation of the data, a neutrino squared 
mass difference in the range: 

Am 2 = 10~ 3 to 1CT 2 eV 2 

and a nearly maximal mixing angle 9 are suggested. The 
interpretation in terms of dominant <-> v T oscillation 
channel is favored (even if it is not possible at present 
to exclude the oscillations into a sterile neutrino). Re- 
sults from earlier experiments on atmospheric neutrino 
anomaly such as Kamiokande ||], 1MB [[|, Frejus 0, 
NUSEX H as well as recent results from MACRO g 
and Soudan-2 are consistent with the results from 
Super-Kamiokande detector. The shape of zenith angle 
distribution of events recorded at Baksan experiment Iq] 
is not in good agreement with the expectation. Also, the 
agreement does not improve much by invoking neutrino 
oscillation scenario. 

Different tests of these observations are essential in the 
perspective of passing from "evidence" to "discovery" . 
Furthermore, new tests should address the problem of 
precise evaluation of the parameters of oscillation. Long 
baselines between the points of neutrino production and 
detection with controlled artificial neutrino beams will 
permit crucial tests, especially if sufficiently large val- 
ues of LIE will be amenable. In this connection, new 
results from Super-Kamiokande, MACRO and Soudan- 
2 @J|,0| will be important to optimize the strategy of 
search (characteristics of the neutrino flux, design of the 
detectors, etc.) and also to interpret the result. At the 
same time, the new neutrino detectors will be able to 
study in more detail the natural (and cost free) atmo- 
spheric neutrinos flux. 

At Super-Kamiokande detector, one of the most im- 
portant observations is certainly the up/down asymme- 
try in ^-induced muon flux. The significance of this 
result is well beyond the level that could be attributed 
to the systematic effects coming from geomagnetic effect, 
detector response etc. The most favored interpretation is 
in terms of flavour oscillations of neutrinos. However, the 
asymmetry has only been observed in f M -induced muons 
with Sub-GeV and Multi-GeV energy range by studying 
events with neutrino interactions inside the detector iden- 
tified as fully or partially contained (FC or PC) events. 
The events of this type are mostly due to low energy neu- 
trino interactions in the GeV range. 

It is very important to extend the search of an up/down 
asymmetry beyond the Multi-GeV region, in order to test 
the interpretation in terms of neutrino oscillations more 
directly, and to further constrain the allowed parameter 
space. This can be accomplished by using the i/^-induced 



muons in the surrounding rock. The advantage of using 
this data sample is the increase in the effective detec- 
tor mass and interaction cross section with the neutrino 
energy, which compensates for the loss in steeply falling 
energy spectrum of neutrinos, roughly as a power law 
(Section |n|). Because of this, the average energy of neu- 
trino for these events is large compared to FC or PC 
events. However, in order to study the up/down asym- 
metry, Ay/m beyond the Multi-GeV energy region, it is 
necessary to obtain the upward and downward going v^- 
induced muon fluxes. In the present study we formulate 
a proposal to achieve this goal by mean of existing data. 

II. UP/DOWN ASYMMETRY 

Many large detectors like Super-Kamiokande and 
MACRO have recorded large numbers of j/^-induced 
muons produced in the surrounding rock and passing (or 
stopping) through the detectors. However, due to shallow 
depth of operation of these detectors, cosmic muon flux 
dominates over the ^-induced muon flux in downward 
direction (except for certain directions in azimuthal an- 
gle, where the mountain is thicker and the shield more ef- 
ficient); hence, detectors cannot distinguish these events 
from each other. The upward ^-induced muon flux can 
be measured accurately, but this is insufficient to study 
the up/down asymmetry mentioned above. However, de- 
tectors operated at deep underground Kolar Gold Fields 
(KGF) B mines can provide the lacking information. 
Due to large depth of operation and flat terrain; be- 
yond a certain zenith angle, the flux of cosmic muons 
is very small as compared to that of ^-induced muons 
in the surrounding rock. Being gaseous detectors, they 
cannot distinguish between the particles moving in the 
upward and downward directions. Therefore, KGF de- 
tectors measure the sum of upward and downward in- 
duced muon fluxes in a given zenithal direction. Hence, 
by combining data from KGF with those from SK and 
MACRO, it is possible to obtain the upward and down- 
ward fluxes. Such an analysis permits the study of the 
up/down asymmetry beyond the Multi-GeV region. 

Some of the essential features of these detectors are 
summarized in Table || It is important to observe that 
the energy thresholds for muons in these detectors are 
different. To get the asymmetry mentioned above, it is 
necessary to have the same cut on visible energy of muon 
in all the detectors. If the energy threshold for KGF de- 
tectors is increased, then there can be a substantial loss 
of statistics. Hence, it is necessary to match the energy 
threshold of other detectors with that of KGF. Inciden- 
tally, this would also lead to further enhancement in the 
statistics of SK and MACRO detectors. 

The extraction of the asymmetry parameter requires 
the results on ^-induced muon flux obtained by the 
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Detector 


Location 


Min. Depth 
(hg/cm 2 ) 


mm 

(GeV) 


Ref. 


Super 


Kamioka 


2700 


1.6 


§ 


Kamiokande 


Japan 








MACRO 


Gran Sasso 
Italy 


3150 


1.0 


I 


Phase- 1 


KGF 


7000 


0.6 






India 






Phase-2 


KGF 


6045 


0.5 


m 



India 



TABLE I. Features of SK, MACRO and KGF detectors. 



SK/MACRO and KGF Collaborations. At present, these 
results are available from the SK experiment for upward 
through going and stopping muons M. MACRO has pub- 
lished the same for through going muons §. The KGF 
experiment has published the event rate of i/^-induced 
muons in rock for one part of the existing data set; but, 
in order to obtain the flux, informations on the angular 
acceptance, and the efficiencies of trigger and detection 
are also necessary. Since the available information is not 
sufficient to extract the proposed asymmetry parameter 
directly from the published experimental data, we have 
evaluated the sensitivity reach of the proposal based on 
certain assumptions. Each of the assumptions mentioned 
below, corresponds to a specific step in the experimental 
analysis. 

• We assume that it is possible to get the j/ M -induced 
muon flux at the same muon energy threshold 
(Efj, > 0.5 GeV), and the same (zenith) angular 
interval from all the detectors. 

• SK and MACRO detectors have recorded a large 
number of upward going i/^-induced muon events 
that are produced in the surrounding rock. More- 
over, SK and MACRO being ongoing experiments, 
will collect more data. Whereas, KGF experiment 
has been stopped in 1992 and is estimated to have 
recorded about 250 events (Section [V) of similar 



kind arriving in upward as well as downward direc- 
tions in a zenith angle cone of 55° < O < 125°. 
Hence, statistical error on downward flux will be 
relatively much higher as compared to that on up- 
ward flux. This is due to the fact that the error on 
downward flux, which is determined by combining 
KGF and SK/MACRO data, mainly depends on 
the statistical strength of the KGF data. Because 
of this, we neglect the errors on the experimentally 
determined upward flux, and treat it as true up- 
ward ^-induced muon flux in our estimate of the 
sensitivity. 

Phase-2 detector at KGF site has observed about 
23,000 cosmic muons jl(J. Similarly, in Phase-1 
detector also a large number of cosmic muons have 
been recorded. These observations are in excellent 



agreement with the predictions of Miyake's empir- 
ical relation jllj up to O < 55° for Phase-1 and 
9 < 60° for Phase-2 detectors. Beyond zenith 
angle of 60° for Phase-1 detector (65° for Phase- 
2 detector) the rate of cosmic muons is negligible 
as compared to the -induced rate. However, in 
the preceding bin, i.e., 55° < 9 < 60° they are 
estimated to be comparable (cosmic muon rate is 
estimated to be about 30% of total rate pig])- In 
order to increase the angular acceptance, we shall 
assume that it is possible to determine the v^- 
induced muon flux from data in this bin, by sub- 
tracting cosmic muon flux obtained using Miyake's 
relation from observed flux in the same bin (and 
similarly for Phase-2 detector in the angular bin of 
60° < 9 < 65° ). 

• Since the aperture area for KGF detector as a func- 
tion of zenith angle is not available, we have as- 
sumed it to be equal for all bins due to nearly cubic 
geometry of the detector. 

• Systematic effects arising due to different detection 
techniques, geomagnetic locations, composition of 
the surrounding rock etc. are not considered in es- 
timating the sensitivity of the proposed method. 

With these assumptions, the sensitivity of the proposal 
to oscillation parameters is obtained by the following 
steps: 



1. 



2. 



The i/^-induced muon flux, ^(Am 2 , 9, is eval- 
uated for each zenith angle bin of 5° in the range 
of 55° < 9 < 125° as a function of two flavour os- 
cillation parameters (for <-> v T oscillation) using 
the available information on neutrino cross section, 
i^-flux, range of muons etc. as described in the 
next section. 

The exposure factor for KGF setup (for preselected 
oscillation parameters) is obtained by normalising 
total flux to the observed number of events as: 



N, 



KGF 



£^ n ^(AmV,90 



(1) 



where 0, = 5° x i+2.5°. The total number of events 
observed at KGF detectors with the suggested cuts 
is estimated to be N^gf ~ 250 (Section IV). 



Using this factor, we calculate the expected number 
of events, NxcFiQi), in each angular bin for the 
KGF (up+down) setup. The same exposure factor 
£ is used to determine Nup(@i), the number of 
upward events in the KGF data set: 



N KGF (Qi) = £x [^(Qi) 
N UP (Q i ) = £x^(Q i ) 



■J>(180° -90], 



(2) 
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To obtain Nyp experimentally, the number of 
SK/MACRO events have to be rescaled to their 
exposure factors, and then multiplied by the KGF 
exposure factor. However, since the statistics of the 
SK/MACRO detectors are quite high as compared 
to that of KGF, the estimated number of upward 
going events can be treated as a true number. The 
sensitivity of the method is controlled then by the 
statistical fluctuations of Nkgf- Using Eqs. (||), 
the total number of events with exclusion of the 
horizontal bin (85° < < 95°), can be obtained. 
We denote them as N KGF and N UP , for KGF and 
SK/MACRO setup respectively. 

4. Now we can evaluate the asymmetry between the 
upward and downward fluxes as: 



av- _ 



N, 



UP 



N 



KGF 



Kp 



(3) 



5. Finally, the sensitivity to the oscillation parameters 
is expressed as the significance of deviation of the 
asymmetry from unity. 

We emphasize once again that the sensitivity has been 
evaluated only on the basis of statistical strength of the 
data. Possible systematic effects, if any, that could arise 
while combining the results from different detectors have 
to be taken into account while performing the compari- 
son. 



III. 



-INDUCED MUON FLUX 



The estimation is based on the current knowledge of 
neutrino flux, its interaction cross section and muon en- 
ergy loss in matter. The ^-induced muon flux in each 
angular bin can be written as: 



d^(Am 2 , 6,E V ,@) = 



(4) 



where we assumed that the ^-induced muon maintains 
the original neutrino direction; an analogue formula holds 
for the P^-induced muon flux. 

The function Y[„ is the muon yield per neutrino. 
The yield increases with energy of neutrino and depends 
on minimum visible energy of muon required by the de- 
tector (E m i n ). The yield is calculated by considering 
the inclusive cross section for muon production dcr^ 
times the number of target nucleons per cm 2 giving rise 
to sufficiently energetic muon above the threshold energy 
Emin (0.5 GeV in the present calculation): 



E » da 
dE, 



-{E^E^x 



M d£ M (5) 



Na is the Avogadro number and R{Efj) is the range of 
muons in the rock. We used neutrino flux from [fl3| . For 
the cross section, we followed the prescriptions of Jl4) , 
splitting (according to the hadronic invariant mass W) 
the quasi-elastic Jig] , the delta-resonance |l6| and the 
deep- inelastic (GRV94 form factors |i7|1 ) contributions 
to the charged current reaction. For the range R{E P ) we 
used the results illustrated in [jr8[ . 

For non-zero neutrino oscillation parameters, the 
flux will get suppressed as: 



dT^{Am 2 ,9,E u ,e) = 

dF° Vll (E u ,Q) x P lu ^(Am 2 ,6,E v ,L u ) 



(6) 



The survival probability P [v _,„ j for <-> v r oscillations 
is given by 

P Wll ^(Am 2 ,e,E v ,L v ) = 

. 2/„«n • i( , „,Am 2 [eV 2 ] ( 7 ) 
l-sm 2 (2*)xsm 2 (o7 L_L_i l j 

where 6 is the flavour mixing angle, Am 2 the square-mass 
difference, L„ is the distance travelled by neutrino be- 
tween its production and interaction point. The distance 
travelled by neutrino, in turn, depends on the zenith an- 
gle: 



(R B + H v y - (i? ffi - H d y sin 2 9 



(R m - H d ) cos 6 



(8) 



i?0 denotes the radius of Earth, H p the height of neutrino 
production in the atmosphere, Hd the depth of operation. 

The estimated ^-induced muon flux .7-^ (Am 2 , 6, O;) 
is obtained in each angular bin by integrating Eq. (^) 
over the allowed range of zenith angle (O; — 2.5° < O < 
Qi + 2.5°) and neutrino energy [E v > E min ). Fig. |l| 
shows the normalised integral flux of ^-induced muons 
as a function of neutrino energy arriving in a nearly hor- 
izontal direction O = 87.5° ± 2.5°. As can be seen from 
the figure; (a) Neutrinos with energies up to 10 GeV are 
estimated to contribute only to 22% of total events (the 
same becomes 12 % at the SK muon energy threshold). 
In presence of oscillations, which mostly affect the lowest 
energy neutrinos, this fraction will diminish further, (b) 
Half of the observed events are expected to be originat- 
ing from neutrinos of energy E v > 50 GeV (for SK muon 
energy threshold, the median is estimated to be at ~ 75 
GeV). This clearly illustrates that the proposed asymme- 
try samples more energetic neutrinos than those in Sub 
and Multi-GeV data sets. It also indicates the necessity 
of comparing data from different detectors at the same 
muon energy threshold as stated in section |l[ 



Na [R{E^) - R(E min )] 
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FIG. 1. Parent neutrino energy distribution for induced 
muons with _E M > 0.5 GeV arriving at nearly horizontal di- 
rection (6 = 87.5° ±2.5°). 



IV. SENSITIVITY OF UP/DOWN ASYMMETRY 
TO NEUTRINO OSCILLATIONS 



The total number of events and asymmetry A^m m U P~ 
ward (95° < 9 < 125°) and downward (55° < 6 < 85°) 
direction are obtained from estimated ^-induced muon 



flux using Eqs. (pJ— 13|). The t^-induced muon flux in up- 
ward and downward direction is expected to be the same, 
except for low energy neutrinos, due to geomagnetic ef- 
fects |]. Hence, the asymmetry will be unity for null hy- 
pothesis, i.e. if neutrinos do not oscillate. Therefore, the 
sensitivity to the oscillation parameters (Am 2 , sin 2 29) is 
expressed in terms of significance of deviation of asym- 
metry from unity, as a function of these parameters. 

The statistical strength of KGF data is the crucial fac- 
tor which controls the significance of a possible deviation 
of the observable A^, /D from unity, and consequently the 
sensitivity to the parameters of neutrino oscillations. The 
number of events recorded by KGF detectors in part of 
their data sample is 213 fllQ| , Using this and the total 
running time of Phase-2 detector || , we estimate the to- 
tal number to be 225 events. The same path length crite- 
ria was applied to Phase-1 and Phase-2 detectors, which 
led to higher muon energy threshold for Phase-1 detec- 
tor (Table |). In order to match the energy threshold 
with that of Phase-2 detector, it is necessary to reduce 
the path length cut by about 20%, which in turn will 
enhance the aperture, and hence the number of events in 



x The SK analysis of East- West effect proves that the theo- 
retical expectations are well met by experimental data 



the Phase-1 detector. In addition to this, there will be 
a further enhancement due to an increase in the angu- 
lar acceptance by an additional 5°. Therefore, in total 
we estimate about Nkgf ~ 250 events from the two 
detectors in KGF site within the zenith angle range of 
55° < 6 < 125°. 

The error on up/down asymmetry is obtained by prop- 
agating the statistical error on KGF observations as: 



^•^U/D — ^U/D 



KGF 



/D N' — N' 

KGF iv C/P 



(9) 



Fig. shows the plot of the asymmetry A^, /D as func- 
tion Am 2 (at maximal mixing). As it can be seen from 
the plot, the asymmetry approaches unity at small or 
large values of Am 2 . At small values, neither upward 
going nor downward going neutrinos get oscillated signif- 
icantly and hence asymmetry will be close to unity due to 
negligible oscillation. For very large values of Am 2 , oscil- 
lation length is quite small as compared to the distance 
travelled by neutrino arriving in the zenith angle range 
of 55° < 6 < 125°. Hence the oscillation probability 
asymptotically approaches to half (at maximal mixing) 
for upward as well as downward directions, making again 
asymmetry closer to unity. This implies that the param- 
eter Au /D cannot distinguish between extreme values of 
Am 2 . 

Average distance travelled by neutrino in upward di- 
rection in zenith angle region of 95° < O < 125° is 
~ 3000 km, and the median energy of neutrino is ~ 50 
GeV (Fig. |IJ). Neutrino oscillations (see Eq. (W) at these 
typical values are maximal for Am 2 « 2 • 10 eV. This 
value is consistent with the best sensitivity point 
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FIG. 2. Up/down Asymmetry beyond Multi-GeV Region, 
assuming maximal mixing angle. 

The significance of a deviation of A^ /D from unity can 
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be obtained from the estimated value of A^ /D , and error 
on it (see Eqs. (§§)). This is translated to a probabil- 
ity, by equating it to the integral of Gaussian probability 
distribution below the estimated significance. We ob- 
tain the sensitivity regions in (Am 2 , sin 2 29) parameters 
space, which corresponds to a chance probability of < 5% 
(and < 1%) to have the estimated value of the asymme- 
try as a result of a statistical fluctuation. These regions 
are shown in Fig. || It can be seen that the proposed 
method is sensitive to Am 2 in the range of 10~ 3 eV 2 to 
1 eV 2 at maximal mixing. This more or less completely 
spans the region of oscillation parameters obtained by 
SK detector to explain their full data sample. 

Sensitivity of the proposed method is derived on the as- 
sumption that experimental data from all the detectors 
will be available at same visible energy threshold. We 
show now that, even releasing this assumption, the con- 
clusions would not change significantly. Using the formal- 
ism described in Section III we estimated that U Pkg f — 
1.22 x UPsk and UPkgf = 1-13 x UPmacro, where 
UPkgf, UPsk and UPmacro denote the upward going 
^-induced muon flux at corresponding energy threshold 
of these detectors in the angular region of interest. If we 
are able to control the error on (UPkgf /UPsk — 1) at 
the level of 17%, or the error on (UPkgf /UPmacro — 1) 
at the level of 28%, then an uncertainty of ~ 3 % would 
be introduced in the expected number of upward going 
events N^ Pl which was assumed to be the true number 
of events in our sensitivity calculations. Since this uncer- 
tainty is small as compared to statistical error on KGF 
data; it would not affect the sensitivity region signifi- 
cantly. It seems to be quite possible to reach the desired 
control on the uncertainty. This is because we only need 
to get the ratio of neutrino induced muon flux at different 
muon energy threshold, wherein most of the uncertain- 
ties on neutrino flux, interaction cross section etc. tends 
to get cancelled out while considering ratios. In this re- 
spect, it is to be noted that ratio of number of stopping 
to through-going muons is predicted to an accuracy of 
13 % at SK 0; the uncertainty being dominated by the 
spectrum of primary cosmic rays. 



V. SUMMARY AND DISCUSSIONS 

Neutrino oscillation signatures are studied in SK de- 
tector by several methods, namely; a) ratio of electron- 
like to muon-like events using contained data sample, b) 
up/down asymmetry of contained muon events, c) shape 
of zenith angle distribution of upward going passing 
through muons and d) the ratio of stopping to through 
going upward muons. It is to be noted that oscillation 
parameters are better constrained by the contained event 
data samples, which corresponds to neutrinos of rela- 
tively low energy. However, these data may not be able 
to improve the constraints further, once the systematic 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Sin 2 (2i?) 

FIG. 3. Sensitivity Region for 2-flavour oscillation 



errors start dominating over statistical errors. Therefore, 
it is necessary to make more effective use of the data sam- 
ple that probes the high energy neutrino spectrum. 

We have shown that it is possible to obtain the 
up/down asymmetry Ay /D , beyond the Multi-GeV en- 
ergy region using the currently available data from KGF, 
SK and MACRO experiments. This is a direct measure- 
ment as it is does not require any a priori knowledge 
of neutrino flux, its zenith angle distribution, interaction 
cross section etc. (the intrinsic asymmetry of the flux due 
to geomagnetic effects only affects the low energy neutri- 
nos, which, however, give a rather small contribution to 
the neutrino parent spectrum shown in Section [II). We 
have demonstrated that this measurement is sensitive to 
the allowed region of neutrino oscillation parameter space 
suggested by the recent results from SK Collaboration [EJ . 
We have not attempted an accurate study of systematics. 
However, we checked that the inclusion of geomagnetic 
effects does not lead to a significant change in the sensi- 
tivity region. 

The best fit value obtained by Super-Kamiokande Col- 
laboration [jlj using Sub and Multi-GeV data is Am 2 = 
3.5 x 10 -3 eV 2 , and maximal mixing; as can be seen from 
Fig. these parameters lie within the sensitive region of 
the proposed method. 



In conclusion, the proposed asymmetry parameter, 
Au /D , entails high energy neutrinos and it could be ob- 
tained using currently available data from different ex- 
periments. The analysis of this parameter permits an 
independent test of neutrino flavour conversion, having 
sensitivity to the range of neutrino oscillation parameters 
suggested by Super-Kamiokande experiments. 
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